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ABSTRACT 


The controlled svstem 1s the ORION satellite spinning about its single axis of svm- 
metry. Hydrazine thrusters are used as the control and are modeled by ideal, constant 
magnitude step functions. 

The system is normalized and driven from non-zero initial angular velocities of the 
two axes other than the spin axis to the final condition of zero. The control profiles 
required to do this are determined based on a desired controller duty cycle. Adaptation 
of the duty cycle changes the ratio of the time the thrusters are on (fuel use) and total 
time to completion of the evolution. 

A comparison between a single axis and a dual axis controller is presented. Simu- 
lation programs for the normalized svstem using a single axis controller, with a 100% 
duty cycle and a varying duty cycle, and a dual axis controller simulation program. with 
each controller having a duty cycle of no more than 50%, are developed. 

The operation of the system is optinuzed using a svstem cost function. An equation 
relating the controller duty cvcle of the dual system to the fuel time trade-off parameter 
of the svstem cost function is required. A nonlinear feedback control algorithm (func- 
tion of attitude angle rates) is developed from iterations of the simulation, and a priori 
knowledge of the form of the control from optimal control theory. This numerical sol- 
ution will allow system designers to incorporate a closed form state feedback control for 


nunimum fuel time strategies using the ORION satellite’s onboard software. 


THESIS DISCLAIMER 


The reader 1s cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er- 
rors, they cannot be considered validated. Any application of these programs without 


additional verification is at the risk of the user. 
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I. INFRODUCTION 


For the last three years a considerable amount of space systems research at the 
Naval Postgraduate School has been devoted to designing a small general purpose sat- 
ellite (ORION). 

The attitude control configuration chosen for the ORION prototype is spin stabili- 
zation. Other attitude control options, three-axis and gravity gradient stabilization, were 
considered for application on the ORION. The pointing accuracy for the gravity gra- 
dient stabilization 1s not highly refined. Although three-axis stabilization can achieve 
high pointing accuracy, it has problems with regard to thermal control, fuel consumption 
and requires a more complex sensor system. [1 : p. 466-479] 

We apply the theory of fuel'time optimized control to a spinning satellite, where a 
single axis of symmetry is assumed. The satellite spins about its axis of symmetry. 
Control, via hydrazine thrusters. maintains or drives the satellite to a zero spin (or an- 
gular velocity) about each of the two axes orthogonal to the axis of symmetry. 

The purpose is to further explore an improved application of spin stabilization to 


the ORION satellite in an effort towards greater fuel efficiency and cost effectiveness. 


Hl. SPIN STABILIZATION THEORY 


A. TWO DIMENSIONAL RIGID BODY 
A simple two dimensional rigid body model serves as the basis for developing the 
necessary equations of motion for the ORION satellite. Figure | portrays such a model 


in which the moment of inertia, / , 1s defined as: 
J = mil? (1) 


where zis the particle mass and / is the inoment arm length. The angular momentum, 


h, relative to the point O 1s expressed as: 


h= Iw =10 (2) 


where @ is the angular velocity and is equal to 0, the rate of change of the angle @. The 


moment of momentum, , of the force, /, about point O 1s: 


M=TFisn@ (3) 


Point O 





Figure 1. Iwo Dimensional Rigid Body Diagram 


This simple two dimensional model clearly states that momentum is a function of 
only three variables. Ilowever, the situation becomes more complicated when a three 


dimensional model 1s considered. 


B. THREE DIMENSIONAL RIGID BODY 

A three dimensional model ts shown in Figure 2 with the three orthogonal compo- 
nents of angular momentum /,, /, and 4, Rotation about any one of the three axes will 
produce a rotation through some angle. For example, rotation through an infinitesimal 
angle, Ad, about the the x-axis will result in two infinitesrmal components of angular 
momentum: hAdk and /,Ag( =) where /, k andi are unit yectors in the Vv, Zz and x 
axes, respectively. Likewise rotating the y-axis through an angle A@ results tn: h,Adi 
and h,A@( — k). And finally, rotating the z-axis through an angle Ay produces: h,Awj 
and h,Aw( — i). The infinitesimal changes in the original angular momentum compo- 


nents, AA, Af, jy and Af,A must also be considered. 


8 (spin) 


o (roll) 





Figure 2. Three Dimensional Rigid Body Diagram 


Recalling that the moment of momentum, or torque, is equivalent to the rate of 
change of the angular momentum, then: 
h 


. A 
Weal AN 4 
aes At (4) 


Adding the components described in the preceding paragraph and taking the limit as Ar 


goes to zero, results in the following momentum equations: 


M,=h,— hy + hé 
M, = hy - hp + hy (5) 
M, =h,—h,O + hyd 


Applving Equation 2 to a three dimensional body gives [2: p. 109]: 


hy Ie ly lez || ¢ 
by |=|—ly ty —hel 6 6 
h, aie —ty, 132 V 


If the products of inertia are ignored [3: p. 51], such that [,=0 for i#j for i and 


j¢ {x,y,z} then Equation 6 reduces to: 


hy, = Lo 
hy = 10 (7) 
h, = 1,. 


The momentum equations can now be written as: 


ML = xh — LyOW + 1,8 
M, =1,,0 — LwWd + 1 .o¢ (S) 
M,= 1.6 — 1,68 + 1,06 


In the case of a rigid body satellite, where external moments are absent [4: p. 524], 


Equation 8 can be simplified by rearranging to become: 


Lx = Uae = 1,,)0W 
byy0 =e Lix)OV (9) 
LVN LO 


If the satellite spins uniformly around the pitch, 6, axis (also called the spin and/or 
y axis) with an angular velocitv of w, and assuming the satellite has a single axis of 
svmmetry, we can conclude that /,,= 1. We will also assume that 0 =w, +e, where ¢, 
the spin error, 1s a small angular velocity error due to perturbation and that ¢ #0, 
b #0, and W #0, with f and W both small. (2: p. 114] 


Incorporating all of these assumptions into the euler equations, they now become: 


ly ad hee . 
= wi 


Hye St (10) 


/,,é = 0 reveals that the rate of change of the spin error is zero and therefore the spin 


error is a constant! Normalize the system by letting: 
i| 
t=(———_) ee (11) 
The state variables are expressed in terms of the system parameters as follows: 
ee o 
=|. (12) 
x} M, 
Substituting Equation 11 into the euler equations and converting to state variable an- 
x, O 1 |] x, 
e\= (1 
X4 —| O X4 
Applving control inputs to the system, the following system equations are obtained: 
8 I cs Om uz, 
‘\s - (14) 
x9 —] Q X> ] e) U> 


Where w, and w, are the normalized thruster control torques for the x, (roll) and x, (yaw) 


notation results in: 


tJ 
ae 


axes, respectively. 
Initially the system will be investigated with uw, set equal to zero and not used. Then 
both uw, and w, will be used at the same time. The following chapter on optimal control 


theory will explore how the values for these two controllers are assigned. 


I. OPTIMAL CONTROL 


Ideal control is a theoretical concept that may not be achievable, since it does not 
account for the physical constraints of the system it is designed for. Conversely, optimal 
control takes the physical constraints of the system into consideration. When a control 
system is chosen based on a given satellite design and a set performance index, opti- 


nuzation is the result. 


A. PERFORMANCE INDICES 

In general, a performance index is a function of system parameters and to a large 
degree defines the character of the optimal control. This, in turn, determines the con- 
figuration of the control system. 

A performance index is normally chosen based upon the svstem’s requirements, 
which mav often be in conflict with each other. For example, one requirement for a 
Satellite system may be to maintain a specified orientation (Within certain given limits 
of accuracv), while at the same time being required to maximize Satellite lifetime. Since 
lifetime is an inverse function of fuel usage, maximizing lifetime means minimizing fuel 
usage. This is in direct conflict with the requirement for maintamung satellite orien- 
tation, which constantly requires fuel for adjustments. The greater the requirement for 
accuracy of satellite orientation, the greater the rate of fuel usage. Consequently, a 
performance index often represents a compromise of system requirements. 

Changing performance indices further complicates matters, since change in per- 
formance indices results in a different optimal control. Practicality also enters the pic- 
ture, since a desired optimal control may be beyond the capabilities of readilv available 
hardware, making the desired optimal control system impractical. 

The difficulty of choosing an appropriate performance index for a complicated sys- 
tem is further compounded by the fact that practicality requires the components of a 
performance index be easily measured or computed through sensors or rate gyros. In 
choosing a performance index, experience has shown that preference should be given to 
the index that is developed from an application rather than one developed from a pure 
mathematical point of view. 

An appropriate control must be chosen which will minimize the performance index, 


J, which is a cost function that is defined as: 


J =f{x(d), to, uD} (15) 


where x(f) is the state vector (Equation 12), z, and 4, are the initial and final umes of the 
systems operation, and u(z) refers to the control vector (Equation 1|4). 

It is important to understand that, in the optimal control application given in this 
paper, J is a function of the control input. u(z) , and the control input is itself a function 
of x. That is, the state vector, x(z}, is included as a parameter for J since x(r) has such 
a great influence upon the control input parameter, u(s). 

The time optimal control system is one of two performance indices to be considered, 


and is defined as: 


yi | 
s=| dt (16) 
f 


where the time interval 7, <1< 14,18 finite. The second performance index to be consid- 


ered is the fuel optimal control svstem: 


ty 2 


T=] > |u(olat (17) 


i=] 


For the purpose of this discussion, u(z) is defined as the thruster’s state at a specified 
time r. Three states will be considered: “off", “on” (positive), and “on” (negative). The 
positive or negative condition for “on” state is an important consideration. as will be 
seen. 

Combining Equations 16 and 17 results in the performance index to be implemented 


in this study. This index is defined as: 


2 
I 

y= fs Ay ufo at (18) 
if) (=| 


where / is a weighting factor that influences the compromise between control response 
and fuel use. 
The weighting factor, /, is a critical parameter that is constant and positive in value. 


By setting 4 equal to zero, Equation 18 becomes a simple minimum time problem. If 4 


is increased to a value approaching infinity, then we have a minimum fuel problem. 
Determining a performance measure for our system reduces to that of determining one 
parameter: 4. The following sections which simulate the normalized system will explore 


this in more detail. 


B. SINGLE VS DUAL CONTROL 
1. One Controller 
a. Minimum Time 

The minimum time switching curve for a spinning satellite is the next topic 
to be considered. When dealing with a control system where only one controller is al- 
lowed to operate (i.e., u, 1s Set equal to zero), u, 1s a function of both states x, and x,. 
When x, versus x, is plotted, regions for the values of u, can be described. This plot 1s 
called the state space or phase plane. In our minimum time problem, where / 1s Zero, 
u, is always turned on. This is the case for an ideal thruster with a non-varving thrust 
magnitude, || = constant. The thruster is always turned on. It can cause the satellite 
body to rotate in a positive or negative direction since the sign of u, can be positive or 
negative. For simplicity, a magnitude of unity for this constant normalized torque will 
be used. The true switching curve for these conditions 1s given in Figure 3 [5: p. 29.]. 
Since this is difficult to model, the switching curve described in Figure 4 1s often used. 
Consequently, the curve in Figure 4 will be used 1n this study when simulating the system 
model when there is only one “on” controller. 

If the minimum time svstem is started at the initial condition: x, = 6.844 and 
x, = —6.844, then the plot given in Figure 5 shows the system’s response. All of the half 
circles above the x-axis have origins at the point x, =—1,x,=0. Likewise all of the half 
circles below the x-axis have origins at the point x, = +1, x, =0. This is not a coinci- 
dence: 

In Figure 4, regressing in time from the origin allows the system to follow 
either of the cusps. The right cusp is chosen (u, = +1), and followed from the origin until 
the x-axis 1s approached at the point x, =+2 and x,=0. In other words, a half circle 
is drawn counterclockwise (reverse time) around the point x, = +1, x, =0 starting from 
the origin to the point x,=+2 and x,=0. When this point is crossed, the switching 
curve boundary is also crossed and the controller changes from u, = +1 to u, =—1, and 
anew arc must be drawn (counterclockwise, as before) starting at the point x, = +2 and 
x, = 0 with a new origin at the point x,=—1 and x,=0. This arc continues until it 


becomes a half or semicircle at the point x, =—4 and x,=0. Continuing in the same 





Figure 3. True Minimum Time Switching Curve for One Controller 





Figure 4. Approximated Minimum Time Switching Curve for One Controller 


ae ; 180 — « 
manner such that two semicircles are drawn above the x-axis and two plus T1830 
G S44 e e ° 
Sonere o = tant ery = 49.506 degrees, semicircles are drawn below the x-axis. By 


simple geometry the final point for this reverse time exercise is the same point as the 





Fieure §. System State Space with One Minimum Time _ Controller: 
Initial Condition (6.844,-6.844) 


mnitial condition for the minimum time (forward time, clockwise) system simulation. The 


conclusion that can be made from this ts that the system is modeled correctly. 


10 


Another check can be emploved to ensure this even further. Since the svs- 
tem is normalized, the time it takes to travel each half circle corresponds to z seconds. 
Therefore in this system model simulation, it is expected that it would take 4.7249 times 
z seconds (14.84 seconds) to go from the inital condition to the zero condition. The 
computer simulation takes 14.82 seconds to get to the point x, = —.QO0000061S. 
x, = —.0005524. Due to the configuration of the computer program, driving the system 
to zero results in iterative reducing of the integration step size of the simulation so that 
much computer time is wasted. This “close enough” final condition strongly implies 
that the system model 1s running as expected. 

b. Niinimum Fuel] Time 

The discussion above dealt with the minimum time problem, where / 1s set 
equal to zero. Now it must be determined how to model the system when the controller 
can be “ofl”, 1.e.. w= 0 and therefore 4 #0. The following paragraphs develop insight 
in the nunimum fuel time control problem using Pontryagin’s Minimum Principle [6] for 
second order linear systems with bounded control inputs. 

Recalling that the svstem in state variable annotation can be written in the 


form of first-order linear differential equations: 
x, =fi{x, u} (19) 


Where X 1s the vector of the state variables (x,,x,), and wis the vector of control variables 
tie). 


The Hamiltonian defined as: 


2 
H = Jug +) phi. W) (20) 
i! 


must be minimized. Here J,,, is the argument (1 + A(|u,| + |u,])) of the cost function, 


J. Also p, and p, are auxilliary variables given by: 


api —6H 

i MES en 
and, 

OS = Web (oR) 


Il 


Substituting Equation 14, Equation 20 becomes: 
H = p, (xo + th) + po( —x, + uy) + 1+ AC la + fel) (235 


Substituting Equation 23 into Equation 21 and solving gives: 


ee: 
a (24) 
| oo mae 
The solution to Equation 24 is of the form: 
= A sin{t — 
Pi (¢— @) (25) 
Pp, = A cos(t — >) 
where 4 and ¢ are integration constants. 
To minimize the Hamiltonian, the form of the optimal control 1s: 
u.=—sign(p,) if lpt>A 
u, =0 if |p| <7 
1 | , hee (26) 
u=-—sign(p,) if |p,l> 
u, = 0 if |p, <2 


The system with no control applied will describe a circle in its phase plane. In other 
words, if the system starts at some initial condition other than zero. the system state 
values will revolve clockwise (forward time) around the origin of the x,, x, state space 
plot. Since no analytical solution was obtained for the minimum fuel’time problem, an 
optimized svstem model is deduced. Based upon Equation 26 and the sinusoidal re- 
quirement described in Equation 25, both w, and u, must cycle from plus to minus every 
180 degrees and be 90 degrees out of phase with each other. Two lines with slopes that 
have the same magnitude but are of opposite sign will achieve the proper periodicity 
requirement, since the regions where the control has a value of zero will always take the 
same period of time to travel through. Figure 6 shows lines with slopes of —1 and +1 
(6 = 90 degrees) and the minimum time switching curve. Since the value for u, above the 


switching curve is —I , and below is +1, there is a conflict in values for the controller in 


2 


the shaded regions A and B. In Table I, the experimental matrix is shown where three 


simulations will be run with varying values for the single controller, 14. 


Table 1. CONTROLLER ASSIGNMENT FOR THE SINGLE CONTROLLER 


Wethod 
Method? [west [w=-1 
Method 





Figure 6. Controller Assignment for One Controller 


Figure 7 shows the curve that 1s common to all three methods. The final port (the one 
closest to the axes’ ongin) for the three simulations is just outside region B, as depicted 


in Figure 6. The separate plots of the three methods from their last common point to 


an end condition described by a circle about the axes’ origin with a radius of .O1 are 
shown in Figure 8. 


The results of the experimental matrix shown in Table | are presented in 
iapley 


Table 2. SINGLE CONTROLLER RESULTS 


Time to End Fuel Consumed 
Condition 

Method 1 10.80 

Method 2 2017 10.75 


Method 3 10.6] 














Ficure 7. One Controller - Minimum Fuel/Time for the Far Field: Initial Condi- 
tion (0.844,-6.844); @ = 90 degrees 


LEGEND 
_ METHOD 1 __ 
METHOD 2 





Figure 8. One Controller - Minimum Fuel/Time for the Near Field: @ = 90 de- 
grees 


2. Two Controllers 
a. Alintnum Time 
The switching curves for the minimum time dual control system are shown 
in Figure 9. The analytical derivation for these curves can be found in [7]. Again, as 


with the single minimum-time controller, these curves are difficult to model, so thev are 
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linearized in the region outside a circle of radius two about the axes’ ongin. This curve 


ieepresented in Figure 10. 





Figure 9. True Minimum Time Switching Curve for Two Controllers 





Figure 10. Approximated Minimum Time Switching Curve for Two Controllers 


b. Ainimum Fuel] Time 

The major difference between the switching curves for the dual control svs- 
tem and the one controller switching curve given in Figure 6 1s that there is now a 
controller on the y-axis. 

The control regions for the case where each of the two controllers has less 
than a 50% duty cycle (Case 1) is depicted in Figure 11. Figure 12 shows the switching 
curves and control areas when the angle @, the deadzone angle, equals 90 degrees. From 
the state equations it can be seen that if x, =0 then 0=x,+u,, or x,=—u,. In other 
words, when x,=+1 then u,.=-—1. Thus the point (0,+1) in the phase plane plot 
equates to the controller value: u,=—1. Likewise, the point (0,—1) corresponds to: 

=+] . This is a subtle but important point. Be aware that the switching curve or- 
gins for the controllers will be located on the x and y axes ONLY when @ = 90 degrees. 
(This will not be the case when @ < 90 degrees. Discussion of this alternative will appear 
later. ) 

Region | in Figure 12 has the following controller assignment u, = —1 and 

= 0; region 2 has vu, = +1 and u, = 0; region 3 has u, = 0 and uw, = +1; and region 4 has 
u,=0 and u,=—1. The minimum time switching curve for the x-axis remains as It was 
for the one controller case; u, has the value of —1 above the curve and +1 below. The 
minimum time switching curve for the y-axis is the new addition. The value for w 1s 
—] to the right of this switching curve, and +1 to the left. 

The four shaded areas in Figure 12 implv a controversial assignment of the 
controller values. Three simulations with varving values for the shaded regions a, b, c, 


and d will be run as described in the experimental matrix shown in Table 3. 


Table 3. CONTROLLER ASSIGNMENT FOR THE DUAL CONTROLLER 
(@ = 90 ecerecs! 


For Se ee 
= = Ja 


Method 
Method 3 






Results from the series of simulations are shown in Figure 13 and 
Figure 14. In Figure 13 the curve that is common to all three methods is shown. The 


final point (the one closest to the axes’ origin) is just outside the region depicted in 


Figure 12 as region a. The plot for the three methods starting from their last common 
point to an end condition described by a circle of radius .Of around the axes’s origin is 


even im Figure 14. Table 4 compares the results of these three methods. 


Table 4. DUAL CONTROLLER RESULTS (6 = 90 degrces) 


Time to End uel Consumed 
Condition 











Method | 10.81 
Method 2 
Method 3 10.60 





Figure 11. Controller Assigument for Two Controllers - Case I: 6 > 90 degrees 


REGION 1 






REGION 4 


wy, 7 x 
b 


Y. | 4 
REGION 2 


Figure 12. System State Space with Two Minimum Fuel/Time Controllers: 
? = 90 degrees 


REGION 3 


20 





Figure 13. Two Controllers - Minimum Fuel/Time for the Far Field: Initial Con- 
dition (6.844,-6.844); @ = 90 degrees 

















N 
\ 
\ 
ON x 
Ps se 
\ 
\ 
\ 
LEGEND | 
METHOD 1 
 MUSTTOD 2 
Oo - 
Xl 
Figure 14. Tyo Controllers - Minimum Fuel/Time for the Near Field: 
@ = 90 degrees 


The different control regions for the dual control case where @ < 90 degrees 
multaneously on. 


(Case 2) 1s shown in Figure 15S. There are four regions where both controllers are s1- 


The switching curve for these minimum time regions 1s shown in 


Figure 10. All other regions comply with the minimum fuel/time switching curves as 
given in Figure 12 using Method 3, from Table 3, for the controller assignment for the 
shaded regions. Case 2 is a much more complicated scenario than previously encount- 
ered. There are now eight regions of controversial controller assignment. This case will 


NOT be used in this study, but 1s mentioned in order to complete the possible scenarios 


for spinning satellite mninimum fuel, time control problems. 





Figure 15. Controller Assignment for Two Controllers - Case 2: 0 < 90 degrees 


3. Analysis 
The switching curves for the single controller with a 100% duty cycle, where y, 


is always on, is shown in Figure 4. This is the minimum time problem. Figure 6 has 


the switching curves for the single controller with the deadzone angle, @, equal to 90 
degrees. The switching curves for the dual controller with a 50% duty cycle for EACH 
of the two controllers, u, and w,. are shown in Figure 12. No more than one of tie 
controllers is ever on at any given time. Table 5 compares the results of these different 


configurations. 


Table 5. SINGLE VS DUAL CONTROLLER RESULTS 


Time to End Fuel Consumed 
Condition 

Single Control: Minimum 14.73 14.73 

Tae =? ¢ =)desrees 

Single Control (Method 3) PANT 10.75 

(= 90 decrees 

Dual Control (Method 3) LOS 10.60 


= 9() decrees 
The end condition for all cases is described by a circle of radius .01 about the axes’ ori- 




















gin. In comparing the niunimum time, single control system to Method 1 of the two 
control minimum fuel-time configuration with @=90 degrees, not listed in Table 5, 
Method 1 achieves 20.6% fuel savings. A fu:2l savings of 28% is achieved when com- 
paring Method 3. Method 3 improves fuel efficiency since there is actually less than a 
50°o duty cycle because the controllers in the four shaded cusp regions (Figure 12 on 
page 20) are turned off. Table 5 confirms that with regard to both time and fuel effi- 
ciency, the dual controller is by far the preferred system. 

What minimum fuel’time configuration of the dual control svstem should be 
employed? This depends on the parameter /, the constant that provides a trade-off be- 
tween fuel usage and time duration of the system being driven from some initial condi- 
tion to a designated end condition. The following chapter develops a possible solution 


to this problem. 


IV NIERICAL ANALYSIS 


The objective is to find a functional relationship between @ , the deadzone angle, and 


A, the fuel’time trade-off parameter. In order to facilitate this Equation I$ is rewritten 


as the following simplified equation: 


J=T+iF (27) 


where 7, the total system time, is 


'f 
T=|"de=y-", (28) 
ly 


and F, the fuel used during the evolution, is defined as 


iy 2 
r= | Sima (29) 


The system runs until the end condition criteria 1s met as defined by 


x? + x5 < 0.01 (30) 


A. NUMERICAL DATA 


G2 


The following steps describe how the data for relating 4 and @ 1s obtained: 


. CHOOSE an initial condition for the dual control program simulation. (The point 


x, = 6.844 and x, = —6.844 is chosen for the first run since this is the initial condi- 
tion used in all of the simulations presented in the previous chapter.) 


RUN the dual controller fuel/time simulation for varving angles of @. (An incre- 
ment of two degrees on the interval from 90 to 180 degrees is used for the first run 
presented below.) OBTAIN the values for 7 and F for each angle of @. 


CHOOSE a value of 7. COMPUTE the values of J for the various deadzone an- 
gles. (A plot of @ versus J for 4 = 15 is shown in Figure 16. Observe the many 
local minima.) RECORD the value of @ and / that gives the minimum computed 
value for J. 


. REPEAT from 3. with a new value of /. 
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Figure 16. The System Cost Function 


Figure 17 shows the first run curve of 2 versus 6. It displays an exponential char- 
acter. Additional values of 4 are iterated using the above steps. Figure 17 is amplified 


with the new data points and the scale is change to semilog. See Figure 18. 
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Figure 17. Linear Plot Relating the Deadzone Angle to the Fuel/Time Tradeoff 
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Figure 18. Semilog Plot Relating the Deadzone Angle to the Fuel/Time Tradeoff 


Parameter 


B. LINEAR REGRESSION 
To determine if the relation between 7 and @ can truly be considered exponential in 
character, a linear regression is accomplished using a commercial computer package 


called Minitab. I-irst, some background is provided. 
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Assume for the moment that the relation between / and @ is exponential in nature. 


Then, this relation can be expressed by the following equation: 
i = Ae” (31) 


Where A and B are both constants. 
Taking the natural logarithm of both sides of Equation 31 results in the following 


new equation: 
In/=InA-+ Be (32) 


A linear regression of the data can now be accomphshed. Minitab computes the 
following constants for the first run simulation data with the initial condition (6.844, 
—6.844): In dA = —5.6465 (or A = 0.00353) and B = 0.070698. 

The goodness of fit is illustrated through the linear correlation coefficient r for the 


number of data points, 7, and is defined as: 


i) In a B. ohn) a In 4 


where 1<i<n. The linear correlation coefficient 1s always between the values of —1 
and +1. Walues of y close to —1 and +1 indicate a strong hnear relationship between the 
variables 6 and In 4, which means that the equation is useful in making predictions for 
a value of A based on a Value of @ or vice versa. 

The value for r based on the data from the initial condition (6.844,—6.844) 1s 0.977. 

Numerical data is obtained for a new initial condition (5,0.001) using the same pro- 
cedure steps outlined above. Ihe values for A, B andr are: A =0.00236, B= 0.073611 
and r= 0.975. 
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V. CONCLUSIONS 


This study provides simulation programs for incorporating an optimized control 
system for a spinning satellite. The developed models support the theorv that there is 
greater fuel efficiency using a dual control rather than a single control configuration. 
Additionally, fuel can be conserved by designing the system response for the maximum 
time permissible for completion of the evolution. 

The software required for this optimal control design is simple, can be easily imple- 
mented and will require very little computer memory. This will allow the ORION to 


operate autonomously while efficiently using the limited onboard fuel reserves. 
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APPENDIX A. PROGRAM FOR ONE MINIMUM TIME CONTROLLER 


This appendix shows the IODE program used for the minimum time single control- 


ler. IODE is an interactive ordinary differential equations package that runs on the 


VM'CMS time sharing system. It was developed at the Naval Postgraduate School by 


Roger R. Hilleary. 


VARIABLES & INITIAL CONDITIONS: 
X1 = 6. 844000000 
X2 = -6. 844000000 


9) 
29) 


2 
a 


SPECIAL FUNCTIONS: 
EAD = X1“*2 + X2**2 


NEAR = ABS(X1) 


Ul mir x2 ,0,1,-1) 


U2 SPE IF(-xIABS(X1)/2+ Xl + X2*ABS(K2)/2,0,1,-1) 


we ODE IF(NEAR,2,U2,U1) 


DERIVATIVES: 
Mei /D(T ) 
X2 


D(X2 /D(T ) 
i U 

D(F /D(T ) == 
ABS(U) 


3} 


;X1 = phi dot 
;42 = psi dot 
;F = fuel 
;T = time 


Send, condition 
<C€riteria 


; absolute value of 
perl 


;if X2 is above the 
;x-axis Ul=-1, else 
gl 


fs pormen( X1,X2) 
; above the switch- 
ine Clive w2—= 1 | 
;else U2=+t1 


;if the absolute 
;value of Xl is 
;less than 2, then 
; U=U2, else U=U1 
 U-= contro! input 


KOMI (oles X2 
tae = =7 1 to 


;F dot = abs(U) 


OUTPUTS: 

TITLE: MINIMUM TIME PROBLEM 

TABULATE. = I X1 X2 F U 
AT INTERVAL . 1O00000000D-01 

PEOR x2 
AGAINST: Xi 

PLOT: Uv GZ U1 
AGAINST: T 

PLOT sen 
AGAINST: T AT INTERVAL . 1000000000 


END CALCULATION WHEN END .LE. .100000D-01 


Sy 


APPENDIX B. PROGRAM FOR MINIMUM FUEL/TIME SINGLE 


CONTROLLER 


This appendix shows the IODE program used for the minimum fuel’time single 


controller. IODE ts an interactive ordinary differential equations package that runs on 


the VM’CMS time sharing system. It was developed at the Naval Postgraduate School 


bw Roger R. Hilleary. 


VARIABLES & INITIAL CONDITIONS: 
X1 = 6. 844000000 
X2 = -6. 844000000 


F= .0 
T= .0 
CONSTANTS: 


P = 1. 000000000 


SPECIAL FUNCTIONS: 
UNEAR = ODEIF( -X1*ABS(X1) /2+X1+X2*ABS(X2)/2,0,1,-1) 


meek = ODEIFCX2+P*X1,0,1,-1) 


UFAR = ODEIF(CABS(X2/X1) ,P,0,-X2/ABS(X2) ) 


oS 


;Xl = phi dot 
;42 = psi dot 
;F = fuel 

;IT = time 

;P = tan(theta) 


; 1f peimt (Xin) 
-15 belowger on the 
; function described 
;by the first 

; argument of the 
ODE TF willie 

; then UNEAR=+1, 
;else UNEAR=-1 


Sl Eee anieaierc:, XZ ) 
;is below or on the 
; line described by 
;the first argument 
ser the OUEIY line 
; then CHECK=+1 
;else CHECK=-1 


‘ite poimes (AL, X2 ) 
;is below or on the 
; lines with slope P 
; AND slope -P 

; OR 

; above lines with 
;slopes P AND -P 
;(lines drawn thru 
fone in) 

; then UFAR=0, else 
; UFAR=-sign( x2) 


UCOND = ODEIF(ABS(X1),2,0,UFAR) ; if magnitude of eal 


;gt 2, UCOND=UFAR, 
; else UCOND=0 


U = ODEIF( ABS( UNEAR+CHECK) ,0,UCOND, UFAR) ,if point (al 


END = X1*X1 + 


DERIVATIVES: 

DI@ GE Bi 616 Ihe) 
X2 

D(X2 /DCT ) 
=x1 +00 

DCF /D(T ) = 
ABS(U) 


OUTPUTS: 


;in REGION A or B 
; then U=UCOND, else 
; U=UFAR 


X2*X2 ;end condition 
;criteria 


X2 


= Xl “doe 
= ;X2 dot = -X1 + U 


;F dot = abs(U) 


TITLE: SINGLE CONTROLLER MINIMUM FUEL? PiMG 


TABULATE: T 


X1 X2 F U 


AT INTERVAL - LOOCCOCUCCD-Ci 


PLOT: X2 


AGAINST: X1 


PLOT: U 


AGAINST: T 


PGOT 


AGAINST: T AT INTERVAL . 1000000000 


END CALCULATION WHEN END .LE. .100000D-01 
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APPENDIX C. PROGRAM FOR MINIMUM FUEL/TIME DUAL 


CONTROLLER 


This appendix shows the IODE program used for the minimum fuel/time dual con- 


troller. IODE is an interactive ordinarv differential equations package that runs on the 


VM/CMS time sharing system. It was developed at the Naval Postgraduate School by 


Roger R. Hilleary. 


VARIABLES & INITIAL CONDITIONS: 
X1 = 6. 844000000 
X2 = -6. 844000000 


F = .0 

t= .0 
CONSTANTS: 

P = 1.000000000 
Q = 1.000000000 


SPECIAL FUNCTIONS: 
MREARL = ODEIF( -X1*ABS(X1) /2+X1+X2*ABS(X2)/2,0,1,-1) 


SaeCKi = ODEIF(X2+P*X1,0,1,-1) 


UFAR1 = ODEIF(CABS(X2/X1),P,0,-X2/ABS(X2) ) 


;X1 = phi dot 

AZ = PS1l dot 

;F = fuel 

;Il = time 

;P = tan(theta) 

;Q = tan(90-theta) 


Tee porntnui Xl yx2) 
;is below or on the 
func rvon 
;described by the 

; first argument of 
;the ODEIF line 

; then UNEAR1=+1 
;else UNEARI=-1 


wikis DOlMe "GNl. v2) 
;is below or on the 
; line described by 
;the first argument 
;of the ODEIF line 
; then CHECK1=+1 
;else CHECK1=-1 


Foabbeemoro pO gheeel( Gil 5.074, 
;is below or on the 
; lines with slopes 
TEeAND =P 

wok 

; above or on lines 
;with slopes P AND 
;-P (lines drawn 
Stnr Onis Poe then 


UCOND1 = ODEIF(ABS(X1),2,0,UFAR1) 


Ul = ODEIF( ABS( UNEAR1+CHECK1) ,0,UCOND1,UFAR1) 


UNEAR2 = ODEIF(X1*ABS(X1)/2-XK2+X2*ABS(X2)/2,0,1,-1) 


CHECK2 = ODEIF(X2-Q*X1,0,-1,1) 


UFAR2 = ODEIF(ABS(X2/X1) ,Q,-X1/ABS(X1) ,0) 


UCOND2 = ODEIF(CABS(X2) ,2,0,UFAR2) 


U2 = ODEIF( ABS(UNEAR2+CHECK2) ,0, UCOND2, UFAR2) 


END = X1*X1 + X2*X2 


DERIVATIVES: 
D(X1 /D(T ) = = 
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; UFAR1=0, else 
; UFAR1=-sign( X2) 


;if magnitude of Al 
;gt 2, UCONDI=UFARI1 
;else UCOND1=0 


;if point ©(x1, we 
;in REGION ager b 
; then U1=0 else 

; V1=UFARI1 


;if point (xXiya 
-is to the leftaen 
;on the function 
;decribed by the 

; first argument os 
; the ODEIF line 

; then UNEAR2=+1 

; else UNEAR2=-1 


;if point (Aly 
;is below or on the 
; line described by 
;the first argument 
;of the ODETR sence 
; then CHECK2=+1 
;else CHECK2=-1 


;if point (X®,2 
;is below or on the 
; lines with slopes 
5Q AND -Q 

Ons 

; above or on lines 
;with slopes Q AND 
;°Q (lines drawn 
;thru origin) then 
; UFAR2=-sign(X1) 
;else UFAR2=0 


;if magnitude of K2 
;gt 2, UCOND2=UFAR2 
;else UCOND2=0 


;if point (Al, xz 
-in KEGION Cc Onea 
; then U2=0 else 

; U2=UFAR2 


;end condition 
;criteria 


;Al dot = X24 


mo + U2 

D(X2 /D(T ) == seenadot = =x)1 + Ul 
-X1 + Ul 

D(F /D(T ) == ;F dot = abs(U1) + 
ABS(U1) + ABS(U2) : abs(U2) 


OUTPUTS: 

TITLE: DUAL CONTROLLER MINIMUM FUEL/TIME 

TABULATE: T X1 X2 Ul U2 F 
AT INTERVAL . 1000000000D-01 

PLOT: X2 
AGAINST: X1 

mei: Ul U2 
AGAINST: T AT INTERVAL . 1000000000 


END CALCULATION WHEN END .LE. .100000D-01 


Ly 
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